Neutrino-induced coherent pion production off nuclei reexamined 
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It is pointed out that so far all theoretical estimates of coherent pion production off nuclei induced 
by neutrinos rely on the "local approximation" well known in photonuclear physics. The effects of 
dropping this approximation are discussed. It is found that in a plane wave approximation for 
the pion, the local approximation overestimates the coherent neutrino-induced pion production on 
nuclei. 
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A precise determination of neutrino oscillation param- 
eters from long baseline experiments requires a precise 
knowledge of the incoming neutrino energy and a de- 
tailed understanding of the neutrino interaction within 
the detection target. Since the targets in all present ex- 
periments involve heavy nuclei, from C to Fe, a precise 
knowledge of the interaction of neutrinos with nuclei is 
obviously required (for a collection of recent references 
see Ref. [l|). In the few GeV region, this involves pion 
production channels, and a detailed study of this inter- 
action is necessary to improve the accuracy of any deter- 
mination of the neutrino energy spectrum. In this energy 
regime, the K2K Collaboration has observed a significant 
deficit in forward lepton direction when comparing their 
measurements with simulations [2|. Since the observed 
deficit is present at forward angles, i.e., small momentum 
transfers, charged current (CC) coherent pion production 
is one candidate interaction that could explain the differ- 
ence. More recent searches, however, found no evidence 
for CC coherent pion production in the low-GeV region 
[3,3- The MiniBooNE Collaboration, on the other hand, 
has reported neutral current (NC) coherent 7r° produc- 
tion [5| at neutrino energies less than 1 GeV. All these 
experimental analyses have in common that fact that the 
coherent fraction is not accessible directly but has to be 
extracted from data assuming specific models for incoher- 
ent pion production; the theoretical models for coherent 
scattering used in the Monte Carlo event simulationsgive 
too large contributions not seen experimentally [!, 0, [H . 

This has triggered a large number of theoretical pub- 
lications on this subject. The theories generally fall into 
two categories. In one class of theories, the partially con- 
serverd axial current (PCAC) is used from the outset, 
and the coherent neutrino-induced pion production is re- 
lated to a pion forward scattering amplitude. This class 
of models assumes that the process is dominated by the 
axial current and that specific nuclear effects play no role 
besides providing nuclear size information. A very recent 
example of this type of approach is the work by Berger 
and Sehgal @. In a fully consistent calculation, the mod- 
els of this first class should emerge as approximations of 
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the second class of models, which start from a theoreti- 
cal description of the nuclear structure and sum the pion 
production amplitude coherently over all target nucleon 
states. A very recent example of this approach is given 
by the work by Amaro et al. Q, which also contains a 
rather complete list of references to both approaches, the 
PCAC based and the ones based on nuclear structure. 

All these studies of neutrino-induced reactions are 
closely related to the experiments on coherent photon- 
induced pion and eta production on nuclei [1, 0, 

MM, 

[l2j ■ While the latter depend only on vector currents, they 
are, of course, closely related to those involving neutri- 
nos in the incoming channel. These reactions were often 
investigated theoretically by invoking the "local approx- 
imation" . This simplifies the hadronic current signifi- 
cantly, because it allows one to pull the A propagator 
out of the pion production amplitude fl3| . which be- 
comes local, and as a consequence allows one to separate 
out the nuclear form factor. Peters et al. [IJ, 0j| have 
discussed the accuracy and the limitations of the local 
approximation for the case of photon-induced coherent 
7r° production as well as for coherent r\ production (see 
also Refs. [H, El, E3 ) • 

The local approximation is the starting point also 
used by the theoretical studies of neutrino-induced co- 
herent nuclear pion production (see Ref. [7| and refer- 
ences therein). In all these studies, the pion production 
amplitude is factorized into a part that contains the pion 
production amplitude and one that contains the nuclear 
size information. Only via this approximation does the 
nuclear form factor emerge in the expressions for coherent 
pion production on nuclei. Since the results of both the 
PCAC-based models and the nuclear physics models rely 
on this factorization, it is of interest to investigate how 
good this approximation actually is for neutrino-induced 
processes. In this brief report, we present the results of 
such an investigation in which we compare the results 
of a calculation using the local approximation with the 
results of a full calculation without this approximation, 
both without pion final state interactions. 

We now briefly discuss the full calculation and then 
show how the local approximation leads to the expres- 
sions commonly used for calculating the coherent pion 
production cross section. 

We assume that the pions are dominantly created via 
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the A(1232) resonance; in Refs. @, [H| it was shown that 
this is a very good approximation for coherent pion pro- 
duction. The hadronic current for a nucleon is then given 
by 
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= iIlc A F(p 2 A )u(pi)k:G a0 (pA)r^(p, q )u(p) ,(1) 

where w-channel diagrams have been neglected. Here /* 
is the ./VA7T coupling constant, F(p 2 A ) a form factor de- 
pending on the invariant mass of the A, k v the pion mo- 
mentum, p' and p the nucleon's final and initial momenta, 
and q the transferred four-momentum. Correspondingly, 
the A momentum is given by pa = p + Q- C A contains 
isospin factors; it is given by 
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where P a p is the usual RS projection operator. The ver- 
tex function in Eq. denotes the standard clcc- 
troweak vertex structure with vector and axial contri- 
butions including the resonance excitation form factors 
for which we take the same set as applied in Ref. [l9| 
(for details see Refs. US El). The propagator [Eq. ©] 
does not include any in-medium changes of the A spec- 
tral function. Their effect is significant, however, they 
affect both, the local approximation and the full calcu- 
lation. In the target nucleus considered here, 12 C, only 
S1/2 and P3/2 states are occupied. The latter reach far- 
ther out than the former and have a higher degeneracy. 
The effect of in-medium changes of the A spectral func- 
tion is thus expected to be smaller in the full calculation 
than in the one using the local approximation where the 
same spectral function is used for all states [cf. Eq. ©]. 

The dyadic product of the currents in Eq. {T]) yields the 
hadronic tensor, which, in turn, determines the v + N — > 
I + N + 7r cross section. 

In coherent pion production, the single particle cur- 
rent [Eq. (pj] has to be summed over all occupied states 
of the target nucleus. This total current then yields 
the hadronic tensor for coherent 7r production on nuclei. 
This is in contrast to incoherent processes, where — in the 
spirit of the impulse approximation — one sums over the 
hadronic tensor of each nucleon. 

Assuming a mean field model for the nucleons gives 
immediately a momentum distribution for the nucleon 
states. This necessitates the appearance of a momentum 
integral in the hadronic current 
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The quantity G a p represents the full Rarita-Schwinger ^j on j g m 
(RS) propagator 



Here the sum runs over all occupied nucleon single- 
particle states in the target nucleus. The single- nucleon 
current on the right-hand side is obtained from that in 
Eq. ([T]) by replacing the free-particle spinors u(p) by the 
momentum representations of the bound-state spinors 
ipi(p) [same for u(j>')\. In the following comparison we 
obtain the latter from a Walecka-t ype mean field model 
for 12 C using parameters from Ref. We stress that in 
this most general expression ((4]) for the nuclear current, 
the momentum integration extends over the A propaga- 
tor as well, since pa = p + q. Because of the presence of 
the A propagator in the single-particle currents and the 
Lorentz structure of the vertex function r^ M , coherent 
pion production does not test the local vector density of 
the nucleus but instead the nonlocal structure of various 
other (tensor) densities. This is the method used success- 
fully to describe coherent photon-induced ir° production 
on nuclei [T3 |. 

The local approximation now consists in fixing the 
momentum of the initial nucleon state in the product 
GaftiPAjT^^ip, q) to some value. An often-used prescrip- 
tion is 
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As a consequence of this "freezing" of the initial nucleon's 
momentum in the transition operator, the momentum of 
the A resonance is also determined. The propagator of 
the A resonance can thus be moved out of the integral 
and then even out of the sum in Eq. This approxi- 
mation basically consists of suppressing the propagation 
of the A resonance and corresponds to the assumption 
of a very heavy A resonance. As a consequence, the 
W, Z + N ^ ir + N vertex becomes local. Thus we are 
left with the current in the local approximation 
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where i runs over all occupied nucleon states and p A = 
(q + k^/2 according to the local approximation ([5]). In 
an r-space representation, this becomes 

(<Z) - 1 & CAF{p ° a) 02 J\- Mr 

m -x p A - m a + iM A r A 

x I ' d^e^-~ k ^ p tr(p(r : r)P aP (p A )T^(p ,q)) . (6) 



Here the trace is taken over the Dirac indices and p(f, r) 
is the diagonal element of the one-body density matrix. 
This is the final result in the local approximation. Equa- 
tion ([5]) shows that the nuclear form factor has been fac- 
torized out; all the other (nonlocal) densities present in 
the full expression no longer appear. 

In the following, we compare the full calculation, based 
on Eq. ([3]) with a propagating A, with the results of 
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FIG. 1: (Color online) CC induced pion angular distribution 
for a neutrino energy of 1 GeV and target 12 C. The dashed 
curve gives the result of the calculation using the local ap- 
proximation [cf. Eq. ©J; the solid curve gives the result of a 
fully dynamic calculation [cf. Eq. Q] . All curves are without 
pion final state interactions. 



the local approximation [Eq. ©] for the target nucleus 
12 C. To isolate the effects of the local approximation, 
both calculations are done in the plane wave approxima- 
tion in which the produced pion is taken to be a free 
particle. Both calculations use the same nuclear struc- 
ture model, i.e., the density and momentum distributions 
are calculated consistently in the same relativistic mean 
field model. For the local approximation, the incoming 
nucleon's momentum has been set to fi = 0, following 
Ref. [lj| . We note that when we use the same approxima- 
tion, we reproduce the local results obtained by Alvarez- 
Ruso et al. [H. 

Figure [1] shows a comparison of the full calculation 
for the angular distribution of the produced pions for 
E u = 1 GeV with the results obtained by using the lo- 
cal approximation. Over a wide angular range there is 
perfect agreement. However, at very forward angles, the 
local approximation gives a cross section that is about 
20% larger than that obtained in the full calculation. 

The difference between the full and the approximate 
calculation is larger at lower energies, as can be seen 
from Fig. [21 where we compare the full and the local 
approximation calculations for the angular distribution 
at the lower beam energy of 500 MeV. Here the difference 
between the two curves is drastic and amounts to a factor 
of w 2 at zero degrees. 

The pion momentum distribution for E v = 1 GeV is 
shown in Fig. [3J It is seen that the local approximation 
overestimates the full result by about 70% at the peak. 
The slight shift downward relative to the fully dynamical 
result is a consequence of the local approximation which, 
as discussed earlier, assumes a very heavy A thus min- 
imizing any recoil effects. We find qualitatively similar 
results for NC induced coherent pion production. 

Recently, Amaro et al. 0] have criticized the assump- 
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FIG. 2: (Color online) CC induced pion angular distribution 
for an incoming neutrino energy of 500 MeV and target 12 C. 
Curves as in Fig. [1] 
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FIG. 3: (Color online) CC induced pion momentum distri- 
bution for an incoming neutrino energy of 1 GeV and target 
12 C. Curves as in Fig. □ 



tion of a vanishing initial nucleon momentum as used 
in Ref. [l9[ (and also in our calculations for the local 
approximation). Therefore, in Fig. [4] we compare the 
results of our full calculations with the results obtained 
by Amaro et al. The latter calculation does not use 
an incoming momentum distribution ~ S(p), but it still 
employs the local approximation from the start. Again, 
the curves shown do not contain any final state interac- 
tions of the pion. As already found above, there is again a 
drastic disagreement between our full calculation and our 
local approximation result (compare the solid and dashed 
curves). For example, in the CC case at 500 MeV (up- 
per part of Fig. [3]), the full calculation is less than 1/2 
the local approximation result at the peak value. The 
result of Amaro et al. (dash-dotted curve), obtained in 
the local approximation, lies significantly lower than our 
"local" result. However, this difference between our local 
approximation and their result is mainly due to the use of 
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FIG. 4: (Color online) Pion momentum distribution for 500 
(CC) and 850 (NC) MeV incoming neutrino energy. The solid 
and the dashed curves are as in the previous figures, the dot- 
ted (double-dashed) ones show the calculation with a modified 
axial coupling for the local approximation (full calculation) as 
detailed in the text; the dash dotted lines are taken from Ref. 
0- 

a significantly reduced axial coupling Cg (0) in the work 
of Amaro et al. [3], as can be seen by comparing their re- 
sult with the dotted line results in which we also reduced 
C£(0) from 1.2 to 0.89 in our calculations. The small 
remaining differences between the two local approxima- 
tion calculations can be attributed to a somewhat differ- 
ent nuclear structure input and the different treatment 
of the incoming nucleon's momentum distribution. The 
double-dashed line in Fig. [4] shows the full calculation 
with the reduced value of C5 . Again, there is a factor 
of rj 2 difference between it and the corresponding local 



curve (dotted). 

So far the plane wave approximation has been used for 
the outgoing pions, while any observable cross section 
contains the effects of the strong pion final state interac- 
tions (FSI) . These are known to not only lower the cross 
section by about 60 % at the lower beam energies 0] , but 
also affect the shape of the pion momentum distributions. 
In coherent events, the quasielastic pion-nucleon scatter- 
ing and the absorption through the A resonance move 
the peak to lower momenta. We thus expect that the 
local approximation result, having more cross section at 
lower momenta, will be somewhat more affected by pion 
FSI than the full calculation result, so that the shapes, 
but not the magnitudes, will become closer to each other 
once pion FSI are included. 

In conclusion, all available calculations for neutrino- 
induced coherent pion production rely on the local ap- 
proximation for the elementary interaction vertex. Only 
this assumption allows one to factorize out the nuclear 
form factor. For the case of photonuclear reactions, 
this assumption had been scrutinized in many theoret- 
ical studies, starting in the 1980s. For neutrino-induced 
reactions, however, we do not know of any previous study 
of this effect. In the present paper, we show that the use 
of the local approximation involves errors in the pion mo- 
mentum distribution that can reach up to 100% and lead 
to different shapes. The discrepancy decreases with the 
neutrino energy (more than a factor of 2 for 500 MeV, 
about 1.7 for 1 GeV, and less for 2 GeV). The differences 
for the pion angular distribution at forward angles also 
decrease with energy. They are of the order of 100% at 
0.5 GeV and 20% at 1 GeV neutrino energy. There is 
thus a general tendency for the local approximation to 
overestimate the coherent neutrino-induced pion produc- 
tion. Full calculations, including pion FSI and medium 
modifications of the A, which go beyond this brief report, 
will yield a final clarification of this point. 
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